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Abstract. Four integrable symplectic maps approximating two Hamiltonian flows 
from the relativistic Toda hierarchy are introduced. They are demostrated to be- 
long to the same hierarchy and to examplify the general scheme for symplectic 
maps on groups equiped with quadratic Poisson brackets. The initial value prob- 
lem for the difference equations is solved in terms of a factorization problem in a 
group. Interpolating Hamiltonian flows are found for all the maps. 



1 Introduction 



Although the subject of integrable symplectic maps received in the recent 

years a considerable attention, the order in this area seems still to lack. Given 
an integrable system of ordinary differential equations with such attributes 
as Lax pair, r-matrix and so on, one would like to construct its difference 
approximation, desirably also with a (discrete-time analog of) Lax pair, r- 
matrix etc. Recent years brought us several successful examples of such a 
construction [1-8], but still not the general rools and recipes, not to say about 
algorithms. 

Recently there appeared for the first time examples where the Lax matrix 
of the discrete-time approximation coincides with the Lax matrix of the 
continuous-time system, so that the discrete-time system belongs to the 
same integrable hierarchy as the underlying continuous-time one (systems of 
Calogero-Moser type [7,8]). We want to present here one more example of 
this type, which can be studied in full (and beautiful) details, - the discrete- 
time analog of the relativistic Toda lattice [9], see also [10-12]. 

The paper is organized as follows. In sect. 2,3 we recall some facts about 
the continuous-time relativistic Toda lattice, its r-matrix structure and the 
solution in terms of a factorization problem in a matrix group. The most 
part of these facts is by now well known, but it turned out to be rather 
difficult or even impossible to find them in literature in the form suitable for 
our present purposes. In sect. 4 we introduce the equations of motion of the 
discrete-time relativistic Toda lattice and discuss their symplectic structure. 
Sect. 5 contains the Lax pair representation for our system, and in sect. 6 
we give the solution of the initial value problem for our system in terms of a 
factorization problem in a matrix group. 

2 Relativistic Toda lattice 

The relativistic Toda lattice with the couphng constant g'^ e M is described 
by the Newtonian equations of motion 



g"^ exp(xjfc+i - Xk) 



■^k-'^k—l 



g'^ exp(xfe - Xk-i) 



1< k< N. 



1 + g'^ exp{xk+i - Xk) 



l + g"^ exp(xfc - Xk-i) ' 



(2.1) 
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with one of the two types of boundary conditions: open-end, 



Xq = OO, Xjv+l = — oo, 

or periodic, 

Xq = Xn, Xn+1 = J'l- 

It is a known fact (although usually not stressed in the literature) that the 
equation (2.1) may be put into the Hamiltonian form in two different ways, 
which lead to two different Hamiltonian functions belonging, remarkably, to 
one and the same integrable hierarchy. 

The first way to introduce the variables pk canonically conjugated to Xk 

is: 

exp(pfc) ij^g2 exp(xfc+i - Xk) ' ^^'^^ 
which leads to the system 

Xk = exp{pk){l+g'^exi){xk+i-Xk)), 

Pk = cxp{xk+i - Xk + Pk) - g^cxp{xk-Xk-i+Pk-i), 

a Hamiltonian system with the Hamiltonian function 

N 

J+ = J2^^P(Pk)i^ + 9^ - Xk))- (2.3) 

fc=l 

The second way to introduce the momenta pk is: 

Xk 

which leads to the system 

Xk = -exp{-pk){l+ g'^exp{xk- Xk^i)), 

Pk = -g'^exp{xk+i-Xk-pk+i)+g^ew{xk-Xk-i-Pk), 

a Hamiltonian system with the Hamiltonian function 

N 

J- = J2 ^^P(~Pk)i'^ + 9^ ^^P(^k - Xk-l))- (2.5) 

k=l 
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The Lax representation and the integrabihty for the flows with the Hamil- 
tonians (2. 3), (2. 5) are dealt with in the following statement. Introduce two 
N by N matrices depending on the phase space coordinates Xk,Pk and (in 
the periodic case) on the additional parameter A: 

N N 

k=l k=l 

N N 

u = J2^kk- >^~^J2 9'^ Xk+Pk)Ek,k+i- 

k=l k=l 

Here Ejk stands for the matrix whose only nonzero entry on the intersection 
of the jth row and the kth column is equal to 1. In the periodic case we 
set £'7v+i,iv = -E'l.AT, = E^y, in the open-end case we set A = 1, and 

En+i,n — En^n+i — 0. Consider also following two matrices: 

T+ = LU-\ T_ = U-^L. (2.6) 

Theorem 1. The flow with the Hamiltonian (2.3) is equivalent to the 
following matrix differential equations: 

L^LB-AL, U^UB- AU, 

which imply also 

t+ = [T+,A], t_ = [T_,i?], 

where 

N N 

A = Z(exp(pjt)+5'^exp(xfe-Xfe_i+pfe_i))£;fcjfc + AZ£;fc+i,jfc, 

k=l k=l 
N N 

B = Y{Qw{Pk) + / exp(xfe+i -Xk+ Pk))Ekk + A X] Ek+i,k- 

k=l k=l 

The flow with the Hamiltonian (2.5) is equivalent to the following matrix 
differential equations: 

L^LD-CL, U^UD- CU, 

which imply also 

f+ = [T+,C], f_ = [T_,D], 
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where 



N 

C = -A"^ Xl5'^exp(xfc+i - Xfc -pfc+1 +pfc)£'fc,fc+i, 

k=l 

N 

D = -X'^^g"^ e^Y>{xk+i- Xk)Ek,k+i- 

k=l 

So we see that either of the matrices T± (they are in fact connected by 
means of a similarity transformation) serves as the Lax matrix for both the 
flows (2.3), (2.5). Note also that the Hamiltonians J± belong to the set of 
invariant functions of T±, as it is easy to check that 

J+ = tr(T±), J_ = tr(T±^). 

It is often convenient to use instead of the canonically conjugated variables 
Xk,Pk another set of variables Ck, dk defined as 

4 = exp(pfc), Ck^ g'^exp{xk+i- Xk+Pk), (2.7) 

which satisfy the Poisson brackets 

{ck, Cfe+i} = — CfeCfe+i, {ck, dk+i} — —Ckdk+i, {ck, dk} — Ckdk 

(only the non- vanishing brackets are written down). In terms of these vari- 
ables the Hamiltonians J± are expressed as 

^ ^ Ck + dk 

J+^^{ck + dk), -^-^I^^:?— ^, 
k=i k=i 

and the corresponding Hamiltonian flows read: 

Cfc = Cfc(4+i + Cfc+i - 4 - Cfc-i), dk = dk{ck - Ck-i) (2.8) 
for the J+ Hamiltonian, and 

Ck — Ck ( J , dk — dk ( ^-J_ j (2.9) 

\dk dk+i) \dkdk+i dk-idk J 

for the J_ Hamiltonian. 
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The fundamental matrices L, U have in new coordinates the form 



N N 

L = ^ dkEkk + Ek+i,k, (2-10) 

k=l k=l 
N N 

U^Y^Ekk- E CkEk,k+i- (2.11) 

k=l k=l 

For the further reference we give here also the expressions in the variables 
Cfc, dk for the matrices involved in the theorem 1: 

N N 

A{c,d,\) = J2idk + Ck-i)Ekk + XjZ^k+iM (2.12) 

fe=l k=l 

N N 

B{c,d,X) = Y.idk + Ck)Ekk + Xj2Ek+i,k, (2.13) 

k=l k=l 

C{c,d,X) = -X~^J2j^Ek,k+i, (2.14) 
k=i "^=+1 

D{c,d,X) = -X-'^-^Ek,k+i. (2.15) 
k=i dk 

3 Algebraic structure 

Here we recall some of the results of [11,12] on the algebraic interpretation of 
the relativistic Toda lattice as a Hamiltonian system on a particular orbit of 
a certain Poisson bracket on a matrix group ([11] deals with a gauge trans- 
formed Lax matrix, which results in a different Poisson bracket on a group). 
The results concerning the difference equations (part c of the Theorem 2 be- 
low) are, to my knowledge, new; however the similar results for less general 
Poisson brackets can be found in [13,14]. 

First of all, we define the relevant algebras, groups and decompositions. 

1) For the open-end case we set g = gl{N). As a linear space, g may be 
represented as a direct sum of two subspaces, which serve also as subalgebras: 
g = g+ © g_, where g+ (g_) is a space of all lower triangular (resp. strictly 
upper triangular) N by N matrices. The corresponding groups are: G = 
GL{N); G+ (G_) is a group of all nondegenerate lower triangular N by 
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matrices (resp. upper triangular N hy N matrices with unities on the 
diagonal) . 

2) For the periodic case g is a certain twisted loop algebra over gl{N): 

g = {r(A) e gl{N)[X,X-'] : Or(A)f]-^ = t{ujX)} , 

where Q — diag(l,a;, . . . ,u;^~^), to — eyiY){2m/N). Again, as a linear space 
g = g+ © g_, where g+ (g_) is a subspace and subalgebra consisting of t(A) 
containing only non-ncgativc (resp. only negative) powers of A. The corre- 
sponding groups are: G, the twisted loop group, i.e. the group of GL{N)~ 
valued functions T(A) of the complex parameter A, regular in CP^\{0, 00} 
and satisfying flT{X)Q~^ — T(a;A); G+ (G_) is the subgroup consisting of 
T(X) regular in the neighbourhood of A = (resp. regular in the neighbour- 
hood of A = 00 and taking the value / in A = 00). 

For both the open-end and periodic cases every r e g admits a unique 
decomposition t — I — u, where / e g+, m G g_. We denote / — 7r+(T), 
u — 7r_(r). Analogously, for the both cases every T e G from some neigh- 
bourhood of the group unity admits a unique factorization T = CU~^, where 
£ e G+, W e G_. We denote the factors as £ = n+(T), W = n_(T). 

Recall also that the derivative dip{T) of the conjugation invariant function 
99 : G I— >■ C is defined by the relation 

, Vii e g. 

£=0 

Theorem 2. 

a) Equip G x G with the quadratic Poisson bracket (38)-(41) from Ref. 
[12], and G with the quadratic Poisson bracket (33) from Ref. [12]. Then the 
set of pairs of matrices {{L{c, d, X),U (c, d, A))} forms a Poisson suhmanifold 
m G X G, the set of matrices {7±(c, d, A)} forms a Poisson suhmanifold in G, 
and the maps (L, U) ^ Tj^ — LU~^ and {L, U) ^ T_ — U~^L are Poisson 
maps from G x G into G. 

b) Let ifi : G ^ C be an invariant function on G. Then the Hamiltonian 
flow on G x G with the Hamiltonian function (p{LU~^) = ip{U~^L) has the 
form 

L = Ln^id^iT.)) - n±{d^{T+))L, 
U = U7i±{dip{T_)) - 7i±{dip{T+))U, 



tr{d(p{T)u) = ^^{Te^'^) 
de 



e=0 
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and the Hamiltonian flow on G with the Hamiltonian function (p{T) has the 
form 

f^[T,n±{dcp{T))], T^T+ or r_. 
These flows admit the following solution in terms of the factorization problem 

e^MT±{0)) ^ c^^t) U^\t), £±{t) e G+, U±{t) e G- 

(this problem has solutions at least for sufficiently small t): 

L(t) = C-\t)L{Q)C_{t) =U^\t)L{Q)U-{t), 

U{t) = C-\t)U{Q)C^{t)^U-\t)U{Q)U--{t), 

so that 

T±{t) = Cl\t)T^{Q,)C^{t) ^Uz\t)T^{m±{t)- 

c) //" / : G H- > G is the derivative of an invariant function on G, then the 
system of difference equations (t & hz) 

Lit + h) = nii(/(r+(t)))L(t)n±(/(T_(t))), 
u{t + h) = ui'{f{T^{t)))um^{f{T4t))) 

defines a Poisson map G x G i— > G x G, and the difference equation 

T{t + h)=Yil\f{T{t)))T{t)n^[f{T{t))), T = T+ or T_ 

defines a Poisson map G G. These difference equations admit following 
solution in terms of the factorization problem 

(r± (0) ) = £± (nh) Uz\nh), £± (nh) e G+ , U± (nh) e G_ 

(this problem has solutions for a given n at least «//(T±(0)) is sufficiently 
close to the group unity I): 

L{nh) = C-\nh)L{fd)C-{nh) = U^\nh)L{0)U-{nh), 

U{nh) ^ C-\nh)U(0)jC4nh) ^U-\nh)U(0)U4nh), 

so that 

T±{nh) = Cl\nh)T±{Q)C±{nh) = U^\nh)T±{Q) U±{nh) . 
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d) The solutions of the difference equations of the part c) are interpolated 
by the flows of the part b) with the Hamiltonian function ^p{T) defined by 

d^{T)^h-Hog{f{T)). 



The part b) of the last theorem explains, in particular, the theorem 1, as 
for J+(r) = tr(r), J_(T) = tr(T-i) we have 

dJ+{T) = T, dJ_{T) = -T-\ 

and it is not hard to check that 

A = 7r+(r_), S = 7r+(r+), C = 7r_(-T-i), D^t:_{-T-^). 

4 Discrete equations of motion and 
symplectic structure 

We propose two systems of difference equations as discrete-time versions of 
the relativistic Toda latti ce. The equations of motion of the first system 
read: 

exp(a;fe(t + h) - Xk{t)) - 1 



exp(a;fc(t) - Xkit - h)) - I 

1 + exp{xk+i{t) - Xk{t))) (l + cxp{xk{t) - Xk-i{t + 
(l + g'^ ex.p{xk+i{t - h) - Xk{t))) (l + g'^ ex.p{xk{t) - Xk-i{t))) 

The equations of motion of the second one read: 

exp(-3;fc(t + h) + Xkit)) - 1 _ 
exp(-a;fe(t) + Xk{t - h)) - 1 

(l + g'^ ex.p(xk+i(t + h)- Xk(t))^ (l + g'^ exp(xfe(i) - Xk-i{t))^ 
(l + 5f2gxp(a;fe+i(t) -Xk{t))^ {l + g'^ e-Kp{xk{t) - Xk-i{t - h 



l<k<N. 
(4.1) 



l<k<N, 
(4.2) 
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The both systems are considered under one of the two types of boundary 
conditions, either open-end: 



Xo{t) = oo, xjv+i(t) = — oo for all t e hZ, 

or periodic: 

Xo{t) = x^it), xi\f^i{t) = xi{t) for all t e hZ. 

The functions Xk{t) in (4.1), (4.2) are supposed to be defined for t G hZ, h > 
(note that the equations themselves do not depend on h explicitly). However, 
it is convenient to consider (4.1), (4.2) as finite-difference approximations to 
(2.1), and in this context Xk{t) are to be considered as the smooth functions 

of i e M. Then the left-hand sides are expanded in powers of /i as 1 + /z— + 

0{h'^), and the right-hand sides as 

l + hl Xk+i——^ J r - Xk-i-—-^ + 0{h ), 

y 1 + g^cxp{xk+i- Xk) 1 + g^ cxp{xk - Xk-i) J 

so that we recover in the continuous limit the equation (2.1). 

Of course, the both systems are closely related, namely by means of the 
time reversion operation. (Note that the underlying continuous-time system 
(2.1) is invariant with respect to this operation). 

We would like to remark that the equations (4.1), (4.2) admit a simple 
non-relativistic limit: set Xk{t) = qk{t) + ct in (4.1) (resp. Xk{t) = qk{t) — ct 
in (4.2)) with c > playing the role of the speed of light; then in the limit 
c — s> oo the both equations tend to one and the same system: 

/ o ^+ h\\ 1 + ^^exp(gfe+i(t) - gfc(t)) ^ , ^ , . 
ex.p{qk{t + h)-2qk{t) + qk{t-h)) = — — ^ -— — -, 1 < k < N, 

1 + g^ exp{qk{t} - qk-i{t)) 

i.e. to the equations of motion of the discrete-time Toda lattice from [3]. 

In the following wc will adopt the notations from [7,8]: if Xk = Xk{t), then 
Xk = Xkit + /i), Xk = Xk{t — h), so that (4.1), (4.2) take the form 

exp(£fc - Xk) - I ^ (1 + g^ exp(a;fc+i - Xk)) (1 + g^ cxp{xk - Xfc-i)) i< k < N- 
exp{xk - Xk)-1 (1 + 5'^ exp(xfc+i - Xk)) (1 + exp{xk - Xk-i)) ' 

(4.3) 
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exp(-Xfc + Xk)-l ^ {l + g^ exp(a;fc+i - Xk)) (1 + g^exp{xk - Xk-i)) i< ^ < jy 
exp(-Xfe + Xk)-l (1 + exp(xfe+i - Xk)) (1 + fi-^ exp(xfe - Xk-i )) ' 

(4.4) 

respectively. Obviously, (4.3), (4.4) are systems of nonlinear algebraic equa- 
tions for Xk, 1 < k < N oi the form 

Fk{x,x,x)^0, l<k<N. (4.5) 

It will be convenient to discuss the solvability of these systems simultaneously 
with the symplectic structure. 

The general recipe to derive the invariant symplectic structure for a dis- 
crete evolutionary equation of the type (4.5) was given in [1,15]: represent 
the equations of motion in the Lagrangian form 

d(^A{x,x)+A{x,x)j /dxk^O, (4.6) 

then the momenta Pk canonically conjugate to x^ are defined as 

Pk^dA{x,x)/dxk. (4.7) 

The map {x,p) i— > {x,p) preserves the standard symplectic 2-form J2k=i dx^/K 
dpk- Note that for the x we have the equation 

Pk ^ -dA{x,x)/dxk, (4.8) 

and then p may be computed as 

Pk^dA{x,x)/dxk. (4.9) 

Remark. Note that if the system (4.5) may be represented in the La- 
grangian form (4.6) with the Lagrangian function A[x, x ), then the system 

Fk{x ,x,x) — 0, 1 < k < N may be represented in the Lagrangian form with 

the Lagrangian function —A{x ,x). Now it is easy to see from (4.7)-(4.9) 
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that the symplectic maps {x,p) {x,p) corresponding to these two systems 
are mutually inverse. This is just the case for (4.1), (4.2), up to minor mod- 
ifications due to the " artificial" parameter h introduced in the definitions of 
momenta below. We prefer, however, to consider all the maps separately, and 
will return to the above-mentioned circumstance at the end of the section 6. 

The Lagrangian functions for the equations (4.3), (4.4) are expressed with 
the help of two functions (f>i{^), 02(0 which are defined by 



0'i(O=log 



exp(0 



h 



0'2(O=log(l + /exp(O). 



It turns out that, just as in the continuous-time case, the Lagrangian func- 
tions and hence the momenta pk may be choosen in two different ways. They 
lead to two pairs of different symplectic maps (one pair for each of (4.1), (4.2)) 
belonging, remarkably, to one and the same integrable hierarchy. Still more 
remarkable, however, is that this hierarchy is the same as in the continuous- 
time case (see sections 5, 6)! 

4.1 System (4.3), the first choice of momenta 

It is easy to see that (4.3) is equivalent to (4.6) with 

N N N 

A{x, X) = J2M^k- Xk) + Y1 M^k - Xk-l) - M^k - Xk-i ). 
fe=l k=l k=l 

Then the definition (4.7) of momenta pk takes the form 

(exp(a;fc - Xk ) -1) 

exp{pk) = r^, (4.10) 

h{l + exp{xk+i - Xk)) 

and the equation (4.8) for x takes the form 

exp(pfc = ^ r-, 4.11 

h[l + 5f^cxp(xfc+i - Xk)) (1 + g-^explxk - Xk-i)) 

(formulas (4.10), (4.11) are to be compared with (2.2)). The equation (4.9) 
for p together with (4.11) implies : 

ex (~ ) = ex ( ) (1 + cxp(.Tfc+i - Xk)) {1 + g^ cxpjxk - Xk^i)) ,^ 
{l + g'^exp{xk+i-Xk)) {l + g'^exp{xk-Xk-i))' 
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To solve (4.3) for x is now equivalent to solving (4.11) for x. It may be 
directly verified that this last equation may be rewritten as: 

/ X 1 +g'^exp{xk+i - Xk) 

hex.p{pk) — ^— = 

1 - exp{xk - Xk) 

1 - exp(xfc_i - Xk-i) 



1 + hexp{pk) + g cxp{xk - Xk-i) 



1 + g'^exp{xk- Xk-i)' 
In terms of the coordinates Ck, dk this means: if we denote 

, . X 1 + - Xk) , . 

ak = hexp{pk) — -. — — , (4.13) 

1 - exp{xk - Xk) 

then 

afc = l + /i4 + l<k<N. (4.14) 

Suppose for a moment that these recurrent relations define Ok as certain 
functions of Cfe, dk- Then, according to (4.13), this means that the equations 
for X are solved. Indeed, it follows from (4.13) that 

/- N Cfc + hck 
exp[Xk - Xk 



hdi. 



To express now the resulting map in terms of the variables c^, dk alone, we 
derive from (4.13) 

1 + g"^ cxp(,Tfc - Xk-i) _ 

— — — ak — riak, 

1 + g-^exp^Xk - Xk-i) 

which together with the previous expression and (4.12) implies: 

^k+i + hck+i J , ak+i — hdk+i /^ir-\ 

Ck = Ck —r , dk = dk — — . (4.15) 

ak + riCk ak — tiak 

Returning to the relations (4.14), we note that in the open-end case Cq = 0, 
hence we obtain from (4.11) the following finite continued fractions expres- 
sions for Ofe's: 

ai — 1 + hdi, 

hci 

a2^1 + hd2+ , . , ; 

1 + ndi 
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OAT = 1 + hdN H 

1 + hdN-1 + 



hCN-2 



l + hdN-2+ hci 

Obviously, we have: 

ak = l + h{dk + Ck-i) + Oih^). (4.16) 

In the periodic case the recurrent relations (4.14) uniquely define the a^'s 
as the iV-periodic infinite continued fractions. It can be proved that these 
continued fractions converge and their values satisfy (4.16). 

Because of (4.16) it is obvious that the map (4.15) is a difference approx- 
imation to the flow (2.8). 

4.2 System (4.3), the second choice of momenta 

The equation (4.3) can be treated also as a Lagrangian equation (4.6) with 

N N N 

A{x, X) ^ - ^(f)i{Xk- Xk) -^(f>2{Xk - Xk-l) + 't>2{Xk - Xk-l)- 
k=l k=l k=l 

Then the definition (4.7) takes the form: 

,n , 2 / (l + /exp(3+i - Xk)) 

""^^""^ = (l-exp(x.-xO) (1+J=exp(..« ^"'"^ 

and the equation (4.8) for x takes the form: 

h{l + exp{xk - Xk-l)) 

exp pfc = ^ , 4.18) 

1 - exp(xfc - Xk) 

(Equations (4.17), (4.18) are to be compared with (2.4)). Equation (4.9) 
together with (4.18) implies:: 

ex (~ ) = ex ( ) (1 +^^cxp(.Tfc+i - Xk)) {1 + g^ cxpjxk - Xk-i)) ,^ 
^ ^ {I + g"^ Q^ixk+i - Xk)) [l + g'^ QW{xk - Xk-l))' 
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As before, to solve (4.3) for x is equivalent to solving (4.18) for x. The 
formula (4.18) may be rewritten as: 

exp(xfe -Xk+ Pk) = exp(pfe) -h-hg^ exp(xfe - Xk-i). 

In terms of the coordinates c^, this means: if we denote 

ffc = g"^ exp(xfe+i - Xfe), (4.20) 

then 

^^dk-h- hVk-i, l<k<N. (4.21) 

Supposing that thes relations define 9fe's as certain functions on Ck,dk, we 
see that (4.20) allows to solve the equation for x. In order to express the 
resulting map in terms of the variables Ck,dk alone, we derive from (4.18), 
(4.20) the relations 

/ ~ X dkdk 1 + (/ exp(,rA;+i - .T-aO hdk 
exp{xk - Xk) - - ' 



Ck 1 + 5-^ exp(xjk+i - Xk) 4+1 
We use them and (4.21) to derive from (4.19): 

Ck + h^k J J dk-hDk-l /Ar,r,\ 

Ck = Ck+i — , dk = dk+i- 7—. (4.22 j 

Returning to (4.21), we immediately obtain in the open-end case Do = 0, 
and the finite continued fraction expressions for Ojt's follow: 



di — 

C2 



d2-h - 

di — n 



_ Cn-1 

d^-i — h 



dN-2 — h — . _ hci 
di — h 
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Obviously, we have 



= ^ + 0{h), l<k<N. (4.23) 
dk 

In the periodic case the recurrent relations (4.21) uniquely define o^'s as the 
A^-periodic infinite continued fractions, that again converge and whose values 
satisfy the relation (4.22). In view of (4.22) it is obvious that the map (4.21) 
is a difference approximation to the flow (2.9). 

4.3 System (4.4), the first choice of momenta 

Turning now to the system (4.2), we see that it is equivalent to (4.6) with 

N N N 

A{x, x) = -Y^ (pii-Xk + Xk) -^hixk- Xk-i ) + Yl 'P'^ixk - Xk-i). 

fc=l k=l k=l 

Then the definition (4.7) of momenta Pk takes the form 

(l-eXp(-Xfe+3)) , 2 / 

/ N {1 + cxp{xk - Xk-i)) 

exp(pfe) = 2 ? ^ , 2 ? nT' ^-24) 

h{l + exp(a;fc+i - Xk)) (1 + g^ exp{xk - Xk-i )) 

and the equation (4.8) for x takes the form 

/ X 1 - exp{-Xk + Xk) , . 

exp pfc) = 2 F ^' (4-25) 

h{l + g^exp{xk+i - Xk)) 

(formulas (4.24), (4.25) are to be compared with (2.2)). The equation (4.9) 
for p together with (4.25) implies : 

expfe) =exp(p>) + - -^'-ji s; - f^'' . (4.26) 

(1 + 5r^exp(xjk - Xk-i)) (1 + 5'^exp(xjt+i - Xk)) 

To solve (4.4) for x is now equivalent to solving (4.25) for x. This last 
equation may be rewritten as: 

/ ~ , X 1 , /X hg^ exp{xk+i - Xk + Pk) 

exp(-Xfe + Xjfc) = 1 -hexp{pk) —z ■ : — ■ 

exp(-a;A;+i + Xk+i) 
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In terms of the coordinates Ck, dk this means: if we denote 

= exp(-5fe + Xfc), (4.27) 

then 

f3k^l-hdk- P^, l<k<N. (4.28) 

Relations (4.28) imply that [3^. as certain functions (continued fractions, 
see below) on Ck,dk- According to (4.27), this means that the equations for x 
are solved. To express now the resulting map in terms of the variables c^, dk 
alone, we derive from (4. 25), (4. 27), and (4.28) the equality 

1 + exp{xk+i - Xk) . hck ^ , , , 
1 + exp(a;jk+i - Xk) Pk+i 

so that (4.26) implies 

/3k — hck-i J J Pk-i + hdk-i r,n\ 

Ck = Ck- 7—, dk = dk———— — . (4.29) 

(3k+i - hck Pk + hdk 

As for the continued fractions expressions for /Sjt's, we have in the open- 
end case: cn = 0, hence 

/^AT = 1 - hdN] 

hCN-l 



Pn-i — 1 — hdN-i 
Pi^l-hdi 



1 — hd^ ' 

hci 



1-hd, ^ 



1 — hdjy 

Obviously, 

Pk^l-h{ck + dk)+0{h^). (4.30) 

In the periodic case the recurrent relations (4.28) uniquely define the /3fe's 
as the A'^-periodic infinite continued fractions. It can be proved that these 
continued fractions converge and their values satisfy (4.30). Because of (4.30) 
it is obvious that the map (4.29) is a difference approximation to the flow 
(2.8). 
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4.4 System (4.4), the second choice of momenta 

The equation (4.4) can be treated also as a Lagrangian equation (4.6) with 

N N N 

A{x, X) = J2 M-^k + Xk) + Yl ^"^i^k - Xk-1 ) - Yl ^2(2;^ - Xk-l ). 
k=l k=l k=l 

Then the definition (4.7) takes the form: 

h{l + exp(xfe - Xk-i )) 

exp(pfe) = ^ — , (4.31) 

exp{-Xk + Xk) - I 

and the equation (4.8) for x takes the form: 

h{l + expjxk - Xk-i)) {1 + g^ exp{xk+i - Xk)) ,^ 
^ (exp(-Xfc + Xk) - 1) (1 + g"^ exp(xfe+i - Xk)) ' 

(Equations (4.32), (4.33) are to be compared with (2.4)). Equation (4.9) 
together with (4.32) imphes:: 

exp(p,) = exp(p,) (l+g'exp(xfe-Xfe_i)) (1 + exp(xfe+i - x^)) ^ 
(1 + g'^ exp(xfc - Xk-i)) (1 + g'^ exp(xjk+i - Xk)) 

As before, to solve (4.4) for x is equivalent to solving (4.32) for x. The 
formula (4.32) may be rewritten as: 

h{l + g'^ cxp(xfc -Xk-i)) _ 
1 - exp{xk - Xk) 

7 , ^ ^ , 2 / , X 1 - exp(xk+i - Xk+i) 

^h + exp{pk)+g cxp{xk+i~ Xk+Pk)— — ^ 7:3 r- 

1 + g^ exp(xjfc+i - Xk) 

In terms of the coordinates Ck, dk this means: if we denote 

1 - oxp(,?7,^i - Xk+i) 

hlk = CkT——^ p r, 4.34 

1 + g^ exp{Xk+i - Xk) 

then 

^ = /i + 4 + /i7fe, 1<A;<A^. (4.35) 

Ik-i 
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This defines 7fe's as certain functions on Ck,dk, wliicli, according to (4.34), 
aiiows to solve tlie equation for x. It is not iiard to derive from (4.34) the 
relations 



(~ \ Ck 1 + 5"^ — j-v-^'Ki-J- -^K) -I , 

exp(a;fc+i - Xk+i) = r-^, — — ^ r = 1 + 



^ _^ fvjk ' l + g"^ exp(xfc+i - Xk) dk ' 

dk 

We use them and (4.33) to express the resulting map in terms of the variables 
Ck, dk alone: 

Ck = Cjfc_i , dk = 4-1:^ ■ (4.36) 

Cfc_i - Ai7fc_i dk-\ + Ai7fc_i 

The expressions for 7fc's in the open-end case follow from (4.35) and 7Ar — 0: 

CJV-I 



dN + 

CN-2 



dN-i + h+ - — — 
d^ + n 



Cl 

7i = 



hC2 

d2 + h^ 



d3 + h+ . hcN-i 
-. + 



dN + h 
Obviously, we have 

7fe = :^ + C'(/i), 1<A;<7V. (4.37) 

dk+l 

In the periodic case the recurrent relations (4.35) uniquely define 7fe's as 
the A'^-periodic infinite continued fractions, that again converge and whose 
values satisfy the relation (4.37). In view of (4.37) it is obvious that the map 
(4.36) is a difference approximation to the flow (2.9). 
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5 Lax representations 



We show in this section that our discrete-time systems admit the discrete 
analogs of Lax representations with the same matrices L,U,T± (see (2.10), 
(2.11), (2.6)) as the continuous-time relativistic Toda lattice. In the following 
theorems we adopt the conventions described before the formula (2.6). The 
dependence of all the matrices below on the discrete time t & hZis supposed 
to appear through the dependence of Ck, dk on t. 

Theorem 3. The symplectic map defined by (4-14), (4-^5) admits a 
representation in the form of matrix equations 

A(t)L(t + h)^ L{t)B{t), A(t)U(t + h)^ U(t)B(t), (5.1) 

so that 

T+{t + h)^ A-\t)T+{t)A{t), T_{t + h)^ B-\t)T4t)B{t) (5.2) 

with the matrices 

N N 

A(c, d,X) = J2 ^kEkk + hXj2 Ek+i,k: (5.3) 

fe=l k=l 
N N 

B{c,d,X)^Y.^kEkk + hXY,Ek+i,k, (5.4) 

k=l k=l 

where a^'s are defined by (4-14), (^nd b^'s by (5.7) or (5.8) below. 

Proof. It is staightforward to check that the matrix equations (5.1) are 
equivalent to the following ones: 

o-kdk = dkbk, hdk + Uk+i = hdk+i + bk, (5.5) 

OfcCfe = Ckbk+i, hck — Ojfc+i = hck+i — bk+i- (5.6) 

Now it would be not hard to check directly that the equations (5.5), 
(5.6) are satisfied with the following identifications: (4.11), (4.12) for dk, dk] 
Ck = 5r^exp(xjk+i - Xk)dk and similarly for Ck] , 

_ _ ^ (1 + exp(xfc+i - Xk)) (1 + exp(xfc - Xk-i)) 

ttk — GXp[Xk ^k) ,9 , \\ /I I 9 / ~ W 

(1 + 5r^exp(xjfc+i - Xk)) (1 + r exp(xfe - Xfe_i)) 
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as in (4.13), and bk = exp{xk — Xk) 

We prefer, however, to work directly with (5.5), (5.6), not using the ex- 
pressions in terms of the variables. Namely, (5.5) is equivalent to 

dk = dk 7-: — , bk = Ofe — — , (5.7) 

ttfc — nak Qk — ndk 

and (5.6) is equivalent to 

flfc+i + hck+i , Ofc+i + hck+\ 

Cfe = Ck — , fafc+i = ak — . (5.8) 

ttfe + hck ak + nck 

The first equations in (5.7), (5.8) coincide with (4.15), and the compatibihty 
of the second ones is equivalent to 

afc(aA-+i - }i(h:+i) _ nfc-i(aA: - /»4) 
Ofe + hck Ofc-i + hck-i 

which is a direct consequence of (4.14). This completes the proof. 

We would like to mention here that the matrices (5.3), (5.4) when com- 
pared with (2.12), (2.13) satisfy 

A(c, d,X)=I + hA{c, d, A) + 0{h'^), B(c, rf. A) = 7 + hB{c, d, A) + 0{h^), 

as it follows from (4.16). 

Theorem 4. The symplectic map defined by (4.. 21), (4-22) admits a 
representation in the form of matrix equations 

L{t + /i)D(t) = C{t)L{t), U{t + /i)D(t) = C{t)U{t), (5.9) 

so that 

T+{t + h) = C{t)T+{t)C-\t), T_{t + h) = B{t)T_{t)B-\t) (5.10) 
with the matrices 

N N 

C(c, d,X) = Y^ Ekk + h\-^ 'kEk,k+u (5.11) 

k=l k=l 
N N 

D(c, d,X)^Yl Ekk + hX-' Yl ^kEk,k+i, (5.12) 

k=l k=l 
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where ^k's are defined by (4-21), and ck's by (5.15) or (5.16) below. 

Proof. It is staightforward to check that the matrix equations (5.9) are 
equivalent to the following ones: 

dk^k = (^kdk+i, dk + h^k-i = dk + hck, (5.13) 

Cfeffe+i = CfeCjfc+i, Cjfc - hVk = Cjfc - hck- (5-14) 

Now it could be checked by means of direct calculation that the equations 
(5.13), (5.14) are satisfied with the following identifications: (4.18), (4.19) for 
dk, dk] Ck = 5'^exp(xfe+i - Xk)dk and similarly for Ck] Ojt = g"^ exp{xk+i - Xk) 
as in (4.2), and 

2 r ~ N (1 ^ exp(,?A^+i - .rA;+i)) (1 + y' exp(,?/,. - .Ta-i)) 



ck = g exp(a;fc+i - Xk)- 



(1 - exp(a:fe + -Xk)) (1 + g'^ exp(xfe+i - Xk)) ' 



Working, however, again directly with (5.13), (5.14), without turning to 
the representation through the Xk variables, we note that (5.13) is equivalent 
to 

T _ , dk — hVk-i _ dk — hVk-i ,^ -I 

"fc — "fc+i-^ r — ) — Ok-; r — ) [o.io) 

and (5.14) is equivalent to 

Ck + hX)k ^ Ck + h^k 

Cfc+i + iiok+i Cfc+i + nok+i 

The first equations in (5.15), (5.16) coincide with (4.22), and the compati- 
bility of the second ones is equivalent to 

^k+i{dk+i — hi)k) _ ^k{dk — hdk-i) 
Cfe+i + hdk+i Ck + hdk 

which is a direct consequence of (4.21). The proof is complete. 

Note that the matrices (5.11), (5.12) when compared with (2.14), (2.15) 
satisfy 

C(c, d,X)^I- hC{c, d, A) + 0{h'^), D(c, d,\)^I- hD{c, d, A) + 0{h'^), 
as it follows from (4.23). 
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Theorem 5. The symplectic map defined by (4-28), (4-29) admits a 
representation in the form of matrix equations 

L{t + h)B{t) = A{t)L{t), U{t + h)B{t) = A{t)U{t), (5.17) 

so that 

T+{t + h)= Ait)T+{t)A-\t), T_{t + h)= B{t)T_{t)B-\t) (5.18) 
with the matrices 

N N 

A{c, d,X) = Y^ akEkk - hXY^ Ek+i,k, (5.19) 

k=l k=l 

N N 

B{c, d,X)^Yl PkEkk -hXY. Ek+i,k, (5.20) 

k=l k=l 

where (3k's are defined by (4-28), and ak's by (5.23) or (5.24) below. 

Proof. It is staightforward to check that the matrix equations (5.17) are 
equivalent to the following ones: 

dk(3k = OLkdk, hdk - Pk-i = hdk-i - (Xk, (5-21) 

Cfc/3fc+i = ttfcCfc, hck + (3k = hck-i + ak. (5.22) 

Now it would be not hard to check directly that the equations (5.5), 
(5.6) are satisfied with the following identifications: (4.25), (4.26) for dk, dk] 
Cfe — fi'^ exp(a;fe+i — Xfe)(ife and similarly for Ck] Pk — exp(— Xfc+Xfc) as in (4.27), 
and 

N (1 + 5'^ exp(5fc+i - Xk)) (1 + g'^ exp{xk - Xk-i)) 
(1 + exp(xfc+i - Xk)) (1 + g^ exp{xk - Xk-i)) 

We prefer, however, to work directly with (5.21), (5.22), not using the 
expressions in terms of the Xk variables. Namely, (5.21) is equivalent to 

T , Pk-1 + hdk-i Pk-i + hdk-i , . 

Pk + hdk Pk + hdk 

and (5.22) is equivalent to 

Pk — hck-i „ /3k — hck-i /r r,^\ 

Cfc = Cjfc- — , ak^Pk+1-3 7—- (5-24) 

Pk+l — n,Ck Pk+1 — It'Ck 
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The first equations in (5.23), (5.24) coincide with (4.29), and the compati- 
bility of the second ones is equivalent to 

Pk+i{Pk + hdk) _ Pk{Pk-i + hdk-i) 
Pk+i - hck Pk - /iCfe_i 

which is a direct consequence of (4.28). This completes the proof. 

We would like to mention here that the matrices (5.19), (5.20) when 
compared with (2.12), (2.13) satisfy 

Aic, d,X)^I- hA{c, d, A) + 0{h^), B{c, d,X)^I- hB{c, d, X) + 0{h'^), 

as it follows from (4.30). 

Theorem 6. The symplectic map defined by (4-35), (4-36) admits a 
representation in the form of matrix equations 

C{t)L{t + h)= L{t)V{t), C{t)U{t + h) = U{t)V{t), (5.25) 

so that 

T+{t + h)= C-\t)T+{t) C{t), T_{t + h)= V-\t)T_{t)V{t) (5.26) 
with the matrices 

N N 

C{c, d,X) = J2 Ekk - hX-' J2 lkEk,k+i, (5.27) 

k=l k=l 
N N 

V{c, d,X) = J2 Ekk - hX-' Yl kEk,k+u (5-28) 
fe=i fe=i 

where 'jk's are defined by (4-35), and Sk's by (5.31) or (5.32) below. 

Proof. It is staightforward to check that the matrix equations (5.25) are 
equivalent to the following ones: 

7fe_idfe = dk - h^k ^ dk- hSk-i, (5.29) 

7fc-iCfc = Ck-i6k, Ck + /i7fc = Ck + h5k. (5.30) 

Now it could be checked by means of direct calculation that the equations 
(5.29), (5.30) are satisfied with the following identifications: (4.32), (4.33) for 
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dk, dk] Ck = g'^exp{xk+i - Xk)dk and similarly for Ck] 6k = g'^ exp{xk+i - Xk), 
and 

_ 2 (~ - ) (exp(-£fc+l + Xk+l) - 1) (1 + (/^ CXp(.Tfc - Xk-l)) 

Ik 9 xp Xk+l Xk ^gxp(-Xfc + Xk) - 1) (1 + g'^ exp(a;fc+i - Xk)) ' 

Working, however, again directly with (5.29), (5.30), without turning to 
the representation through the Xk variables, we note that (5.29) is equivalent 
to 

T , dk + h^k . dk + h-fk rro-,\ 
dk^dk^i- — , ()jt_i = 7fe_i- — , {b.6l) 

dk-i + h^k-i dk-i + h^k-i 

and (5.30) is equivalent to 

Ck - h^k . Ck - h^k oo^ 

Ck-i - n'jk-i Ck-i - li-fk-i 

The first equations in (5.31), (5.32) coincide with (4.36), and the compati- 
bility of the second ones is equivalent to 

7fc(4+i + ^7fc+i) ^ 7fc-i(4 + fe7fc) 
Ck - h'jk Ck-i - h'jk-i 

which is a direct consequence of (4.35). The proof is complete. 

Note that the matrices (5.27), (5.28) when compared with (2.14), (2.15) 
satisfy 

C(c, d,X) = I + hC{c, d, A) + 0{h^), V{c, d,X) = I + hD{c, d, A) + 0{h'^), 
as it follows from (4.37). 



6 Factorization problems and interpolating Hamil- 
tonians 

It is very remarkable that the matrices A, B, C, D and A, B, C, V from the 
previous section may be identified with the certain factors n-|-(/(r-t)), as in 
the Theorem 2. 
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Theorem 7. There hold following relations: 

A{c,d,X) = U+{I + hT+{c,d,X)), (6.1) 

B(c,d,X) = U+{I + hT_(c,d,X)), (6.2) 

C{c,d,X) = UZ^ (l -hT-\c,d,X)) , (6.3) 

D(c,d,A) = UZ^ [l -hTz\c,d,X)) . (6.4) 
Proof. Define following two matrices: 



Q-{c,d,X) = ^Ekk-X ^ ^ — Efe^fc+i e G_, 
fc=i fe=i 

(5+(c, A) = ^ — Ekk + Ek+i,k € G+- 
fe=i ^'i' ifc=i 

Note now that the recurrent relation (4.14) is just equivalent to the matrix 
equality 

C/(c, d, A) + hL{c, d, A) = A(c, d, A)g_(c, d, A), (6.5) 
and the recurrent relation (4.21) is just equivalent to the matrix equality 

L(c, d, A) - hU{c, d, A) = Q+(c, d, A)D(c, A). (6.6) 

Multiplying (6.5) from the right by U~^, we obtain: 

I + hT+^ AQ_U-\ 

and since Q-U~^ e G_, we obtain (6.1). Prom the previous equation with 
the help of (5.1) we derive also 

/ + /iT_ = U-^AQ_ = BU-^Q_, 

which proves (6.2) in view of U^^Q^ G G_. 

Next, multiplying (6.6) from the left by L~^, we obtain: 

I - hTZ^ = L-^g+D, 

which just implies (6.4) because of L~^Q^ e G+. Finally, from the previous 
equation we derive with the help of (5.9): 

/ - hT~^ = g+DL-i = g+Z-^c, 
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which means the vahdity of (6.3) because of Q^L~^ e G+. 
The theorem is proved. 

Theorem 8. There hold following relations: 

A{c,d,X) = n;^ ((/-/iT+(c,d,A))-^) , (6.7) 

B{c,d,X) = ((7-/iT_(c,d,A))-^) , (6.8) 

C(c,d,A) = n_((7 + /ir-i(c,d,A))-^), (6.9) 

V{c,d,X) = U_(^{I + hTz\c,d,X))-^) . (6.10) 

Proof. Define following two matrices: 

P.{c,d,X) = Y.Ekk-X-'Y.J^Ek,k+ieG_, 

k=l k=l f^^^+l 

JV N 

P+{c,d,X) = J2^Ekk + Xj2Ek+,,keG+. 

k=l Ik-l 

Note now that the recurrent relation (4.28) is just equivalent to the matrix 
equality 

f/(c, d, A) - hL{c, d, A) = P_(c, d, X)B{c, d, A), (6.11) 
and the recurrent relation (4.35) is just equivalent to the matrix equality 

L(c, d, A) + hU(c, d, A) = C(c, d, X)P+(c, d, A). (6.12) 

Multiplying (6.11) from the left by U~^, we obtain: 

and since U~^P- e G_, we obtain (6.8). Prom the previous equation with 
the help of (5.17) we derive also 

I-hT+ = P_BU-^ = P_U-^A, 

which proves (6.7) in view of P^U'^ e G_. 

Next, multiplying (6.12) from the right by L"^ , we obtain: 

/ + hT7^ = CP+L-\ 
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which just imphes (6.9) because of P+L ^ G G_|_. Finally, from the previous 
equation we derive with the help of (5.25): 

I + hTZ^ = L-^CP+ = VL-^P+, 

which means the validity of (6.10) because of L^^P+ G G+. 
The theorem is proved. 

Substitute now the expressions (6.1)-(6.4) into the discrete Lax equations 
(5.1), (5.2), (5.9), (5.10), and (6.7)-(6.10) into (5.17), (5.18), (5.25), (5.26). 
One recognizes immediately the difference equations from the part c) of the 
Theorem 2. This gives us the solution of the initial v alue problem for the 
dynamical system (4.14), (4.15) in terms of the factorization of the matrices 

{{I + hT^mr , 

the solution of the initial value problem for the dynamical system (4.21), 
(4.22) in terms of the factorization of the matrices 

{{I - hT^\o))y . 

, the solution of the initial value problem for the dynamical system (4.28), 
(4.29) in terms of the factorization of the matrices 

{{i-hT^{o)rf. 

, and the solution of the initial value problem for the dynamical system (4.35), 
(4.36) in terms of the factorization of the matrices 

{{I + hT^\o))-'y . 

The part d) of the Theorem 2 may be formulated in our case in the 
following lines. 

Corollciry. The interpolating Hamiltonian for the map ((4- 14), (4-15) is 
given by 

<p+{T) = tr($+(T)) = J+(r)+0(/i), where $+(0 = h-' ^ log(l+H, 

Jo T) 
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for the map (4. 21), (4.22) - hy 



99_(T)=tr($_(T)) = J_(T)+0(/i), 



where $-(0 



h 



-1 



°° d-q 



log( 



1 



1 — hr) 



for the map (4.28), (4.29) - by 

MT) = tr(*+(r)) = J+(r)+0(/i), where *+(0 = h'' f ^ log - 



We have introduced the difference approximations for two Hamiltonian flows 
from the relativistic Toda hierarchy. They turned out to belong to the same 
hierarchy. The inclusion in the general scheme of symplectic maps on groups 
equiped with quadratic Poisson brackets allowed to solve the difference equa- 
tions in terms of factorization problem in the group and to find the interpo- 
lating Hamiltonians. 



and for the map (4-35), (4-36) - by 




where 




7 Conclusion 
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